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Abstract 

The K — > tttt system is analyzed in the chiral limit within the Standard Model. 
We discuss how to connect the short-distance running in the |AS| = 1 case to the 
matrix-elements calculated in a low-energy approximation in a scheme-independent 
fashion. We calculate this correction and the resulting Wilson Coefficients. The 
matrix elements are calculated to next-to- leading order in the 1/N C expansion and 
combined with the Wilson coefficients to calculate the two isospin amplitudes and e' K . 
The AI = 1/2 rule is reproduced within expected errors and we obtain a substantially 
larger value for e' k /ek than most other analysises. We discuss the reasons for this 
difference. We also suggest that the X-boson method is an option for lattice QCD 
calculations. 
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1 Introduction 



The study of CP-violation took a major step forward in 1999 with the first direct evidence 
by two independent experiments for direct CP- violation JTJ corroborating the earlier NA31 
result. The average of the two experiments 
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(2.13 ±0.46) • 10" 



(1) 



is higher than the standard calculations obtained p], || and this has caused claims for 
physics beyond the Standard Model, see e.g. |3J. The purpose of this paper is to calculate 
this quantity within our earlier proposed method and see whether the above number can 
be reached or not. Other opinions critical of this claim are expressed in |J. 

CP-violation was first observed in the sixties and was then shown to proceed only via 
K°-K° mixing. This is called indirect CP-violation and is described by the parameter ek 



A[K L -> (7T7r) J= o] 
A[K S - (7T7r) /= o] 



In addition there is also direct CP-violation in the decay-amplitude possible in K 
decays. This is parameterized by 
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A[K L - (7T7r) z=2 ] A[K L - (7T7r) /=0 ] A[K S - (7r7r) J=2 ] 



A[K S - (7T7r) /=0 ] A[K S - {w) I=G \ A[K S 
Combining both definitions the ratio is directly accessible from 



(7T7r)/=o] 



(2) 



(3) 



A[K L -> (7T7T) J=2 ] A[K S -> (7T7T) J=2 ] 

A[K L ^ (7T7r) /=0 ] A[K S (™) I=0 ] 



(4) 



This last formula is the basis for the double-ratio method used in the actual experiments 

Theoretical calculations of £' k /£k basically reduce to calculating the K — >• itn ampli- 
tudes. This is in the Standard Model a several step process, see Section | and e.g. the 
lectures of for more details. First the heavy particles, top, Z and W, are integrated out, 
resulting in an effective action in terms of four-quark operators. Then the renormalization 
group evolution equations are used to bring the coefficients in this effective action down 
to a hadronic scale thus resumming large logarithms containing heavy masses. This part 
has been performed including gluonic and electroweak Penguin-diagrams to one-loop in 
and to two-loops in ||. 

This has then been used with matrix-elements from the lattice and from the first 1/N C - 
methods M to calculate £' k /£k in @- We will not comment further on the lattice calcu- 
lations. A recent suggestion about their calculation and more references can be found in 
T0| . We will concentrate on the 1/N C method of calculating the matrix-elements. 
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The l/A c -method was first proposed in |9j after it was noticed that using large N c 
instead of vacuum-insertion-approximation for the matrix elements gave a significant im- 
provement for the A J = 1/2 rule. They simply identified the cut-off in meson- loops with 
the scale in the renormalization group. A better way to provide this identification was 
given in [IT, [12|] by using colour- singlet bosons to provide the identification of scales. To 
one-loop it was proven in |12] that this procedure gives the correct matching. The work of 
|13| is essentially the continuation of |§ using this identification directly with the output 
of the renormalization group. 

Our work has concentrated on two more improvements. We realized that our method 
of using colour- singlet bosons, hereafter referred to as the X-boson method, allowed both 
a solution to the scale-identification and to the scheme- dependence that appears at two- 
loops in the renormalization group. This was argued in [IJJ] where we used available results 
to include part of the needed correction. We proved this explicitly in the case of Bk in 



and below for the IAS* 



1| case. The other improvement we have included is the 
use of the ENJL model for the couplings of the X-bosons to improve on the high-energy 
behaviour. For Bk this can be done directly and we presented results including all quark- 
mass corrections in |12|, 15| . For K — > titt decays this produced too unwieldy calculations 
and we use instead the method of [|16| to reduce the calculation to off-shell two-point 
functions. Results for the operators Q\ to Q§ in the chiral limit were presented in |L4] and 
in the same limit for Qj to Qio below. 

Related work, concentrating on a more QCD inspired improvement of the X-couplings 
is Ref. 
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We regard this method as very promising but it has not yet reached the stage 
where results for the quantities considered here are possible. 

An approach where no attempt at precisely identifying the scales in the matrix-element 
and the short- distance running is done is the chiral quark model approach of On the 
other hand this series of papers has attempted to include several effects neglected by the 
standard analysises || by systematically going to p 4 and thus provided an indication that 
the errors in |2| were severely underestimated. 

Another approach to determining the non-leptonic matrix-elements is to use QCD sum 
rules to calculate more inclusive quantities and then relate those to the matrix-elements. 
Examples of this approach are [19[ . 



As discussed in [[12| for the case of the Bk parameter, a general problem with this 
method is the existence of operators in the CHPT description that contribute to the inclu- 
sive quantities but not to the decays K — > tttt. The reason is that the so-called spurion-fields 
have now non- vanishing derivatives Q. 

We would like to stress that our method basically reproduces the AI =1/2 rule, has 
no free input parameters as several other models have and includes a correct scheme and 
scale identification at all stages. As will become obvious, we also reproduce the value of 
£ 'k/ £ k within the uncertainty of our approach. 

Our X-boson method can be used in general and is not restricted to 1/N C . E.g., 
by adapting the method of |2TJ to use A-bosons instead of photons it might provide an 

A 33 )). 



2 In the notation of Section || an example of an extra octet term is tr(A 32 u AI )tr(V A1 (An + A 2 
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alternative lattice QCD method to calculate K — » im on the lattice. 

In the next section we introduce the basic notation used to describe K — > txtx decays 
to lowest order in Chiral Perturbation Theory (CHPT). Section |3] describes the X-boson 
method and calculates the corrections needed to identify the scales in QCD and the low- 
energy model used. The main result is presented in the appendix. Section |] presents the 
calculation in 1/N C using CHPT and the ENJL model for the X-boson couplings for the 
electroweak Penguin operators. The numerical input used is summarized in Section [5[ The 
results from the previous sections and the matrix-elements of [H| are then combined to 
calculate the parameters in the CHPT AS = 1 Lagrangian in Section ^ and presented 
using the more standard bag parameters, Bi in Section [7|. The next section contains the 
calculation of e' K and briefly of Ek as well. The last section discusses the result for e' K /eK 
and adds the two corrections which are expected to dominate. 



2 Chiral Perturbation Theory Lagrangian 

We use the standard Chiral Perturbation Theory (CHPT) counting with p a momentum, 
energy or meson-mass and e the electromagnetic coupling constant e = \fAixo>. The lowest 
order strong and electromagnetic Lagrangian is given by 



£(2) = _JL tr fatf + x+) + e 2£ 2tr (QUQrf 



(5) 



with U = u 2 = exp(A a 7r a / F Q ) and = iu^D^Uu^ . X a are the Gell-Mann matrices and the 
7i a are the pseudoscalar-mesons ir, K, and 77. Q — diag(2/3, —1/3, —1/3) is the light-quark- 
charge matrix and x+ = u^x u ^ + ux^u. x = 2-Bo diag(m u , ma, m s ) collects the light-quark 
masses. To this order F = F n is the pion decay coupling constant. An introduction to 



CHPT can be found in the lectures pl |. 

To order e°p 2 and e 2 p°, the chiral Lagrangian describing \AS\ = 1 transitions is given 

by 



,(2) 

'|AS|=1 



CF VG £ tr A 32 g + CF { 



G 27 f^tr(A 4J ^)tr(A H ^) 



G 8 tr {A^uX) + G" 8 tr (A 32 ^ 



h.c. 



with A i:j = u\ i:j u\ (Xij)ab = $ia8jb, Q = u ] Qu 
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v ud v: 



-1.06- 10~ 6 GeV" 
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The SU(3) x SU(3) tensor t ij < H can be found in 



The K — > 7T7T invariant amplitudes can be decomposed into definite isospin quantum 
numbers amplitudes as [A = —iT] 



A[K S 
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Where K s ~ K[ + e K K$\ K° 1{2) = (K° - (+)iT)/V2, and CP(K 1 ° (2) ) = +(-)K 1 ° (2) . The 
final state interaction phases 5q and 5 2 are included into the amplitudes Aq and A2 as 
follows. 



A, = 
A 2 = 



-1 ao e 
-i a 2 e 



(do 



(52 



for the isospin 1/2 amplitude; 
for the isospin 3/2 amplitude . 



(9) 



Using the Lagrangian of Eq. (pj) we obtain 



a 



(1-2 
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CF n 



10G 27 (m 2 K - m^) - 6e 2 G* B F ( 



(10) 



In the presence of CP- violation the couplings Gg, G27, and Ge get an imaginary part. In 
the Standard Model, ImG27 vanishes and ImGg and ImG^ are proportional to Imr with 
r = -Xt/Xu and A* = V^V^. 



3 Short-Distance 

In Section f| we will calculate the long-distance contributions in the X-boson scheme. This 
scheme requires only the matching of currents between the low-energy model/theory and 
the short-distance QCD calculations. This way the scheme- and scale-dependence inherent 
in most of the other calculations is consistently avoided. This scheme has been shown to 
one-loop to reproduce the correct matching [1T2J and was argued to work to all orders in 



||14| . The precise proof to two- loops was done in |TJ] for the case of the 5^-parameter. 
Exactly the same calculation can be performed for the \AS = 1| case. There is no new 
problem but in practice the calculation is much longer because there are now ten operators 
instead of one. Penguin diagrams need to be taken into account and photonic loops as well 
as gluonic ones have to be considered. We will therefore present fewer explicit expressions 



than in 15 



We start from the effective action that is the output of the renormalization group 
running at two-loops which is given by 

G 
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with Cj — Zi + yiT and 
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with {qi^q')(L,R) = 97^(1 T 75)9'; a an d P are colour indices. 

The Wilson coefficients Cj(//r) are dependent on the scheme chosen for 75, the choice 
of evanescent operators, the scale in the renormalization group as well as other scheme 
choices like the choice of infrared regulators. All of them have to be consistently treated 
when calculating matrix elements in order to reach a physical result. Reviews of this 
can be found in || and 0. The main original references are ||. We have chosen two 
slightly different implementations of the running, both in the NDR scheme. One where 
we linearize fully the NLO part, this way all divergences appearing in the expressions 
can be analytically resolved 3 and the scheme dependence is canceled exactly except for 
effects of order as(m c ) 2 . The other implementation is where we exactly solve the two- loop 
evolution equations for the Cj(/xr). In this case the scheme-dependence is only canceled 
to order as(f^R) 2 . We regard the difference between these two as a rough estimate of the 
uncertainty due to higher order corrections in ois(//r). 

The next step is now to replace the effective action flTT| ) by another equivalent effective 
action in D = 4 where the quarks only appear in vector, axial-vector currents, and scalar 
or pseudoscalar densities. These we know how to hadronize. Specifically we replace Eq. 
(O) by 

Tx = 9i x i ((«7/A + (w7m«)l) + 5-2^2 ((s7/^)l + (^tA) 

+#3^3 I {sind) L + E (97M9)i)+^4 E X qA((s-f^q) L + (w^l) 

\ q=u,d,s J q=u,d,s 

3 Note that in none of the papers of Q all necessary formulae were given, we have worked them out 
ourselves. 
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q=u,d,s 



-2){qd) 
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q=u,d,s 



+9rX¥ + I] o e q(Q^Q)R ) + 9s E X g,& \( s i)l + (-^)7,e q (qd) R 



q=u,d,s 



q=u,d,s 



+g 9 Xg I (s7 M d) z + 51 7^s(97^)l J + #io E -X£io + -e q {qi^q)i\ . 

\ q=u,d,s J q=u,d,s 

(13) 

Here all colour sums are performed implicitly inside the brackets and (qq')(L,R) = ?(1 = F75)? / - 
For simplicity we choose all X-bosons to have the same mass. We now determine the 
couplings Qi as a function of the Cj by taking matrix elements of both sides between quark 
and and gluon external states, labelled by ip and if)' and require 
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Both the left- and right-hand-side can be written in terms of the tree level matrix-elements 
of the Qi between the states if), if)'. Leading to 



+ 



E 

3 

a | 

7T 



S i:j + 



a s (fi 



R 



m 



7T 

m 



7ij log + ry + F^m 2 , if), if)') 



= £ 
- ^(( 



5y + 



+ r?. + F?.(m 2 ,^^) 



7T 



(7»j - 7y) lo 



m „ m _, , o .A 

g tt + To" log — + r„ + Fy(m , ip, ip ) 

M x fie ) 



Ti 



7ii 



log 



m 



m 



+ 7" log— +r! j +F°Jm 2 ,if),if)') 



M x ' 113 to nc 13 13 



9j{^ 



c 



Ml 



(15) 



The parameter m is some IR regulator. For conserved currents 7^ = and there are 
no logarithms of the scale \xc- For non-conserved currents or densities 7^ can be related 
directly to their anomalous dimensions. 

The infrared sensitive parts are in the functions F^ and in Ff- and are identical on 
both sides. We have checked this explicitly for a set of states similar to the one chosen in 



15 1 . The resulting formulas are very long and only reinforce the argument for complete 
infrared ambiguity cancellation. Notice that the matrices 7^ and 7?- are one-loop anoma- 
lous dimension matrices of the set of four-quark operators ( |12D [0. The left-hand-side is 
thus completely scheme-independent to order a 2 s . The scheme-dependence of the Cj(fj,ji) 
is canceled by the scheme-dependence of the constants and rfy The right-hand-side is 
similarly scale-independent 

From Eq. (|T5|) we can now obtain 
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a ( M x _ P . Mc P ~ P \ 

M 7&log 7^ +7 « log ^ +r *- r «) 

+0 (a s (fi R ) - as(^c))] Cj(/*Ji) • (16) 

The X-boson couplings are now completely scheme- and scale-independent and the large 
logarithms of My/ are resummed and included in the values of the couplings g\ . The explicit 
expressions for the differences — fy and rfj — rfj are given in the appendix. 

We have not used any large N c arguments up to this point. Now we need to calculate 
the matrix-elements of X-boson exchange between hadronic external states. As described 



in more detail in |Lo| and references therein, we do this by rotating the integral over the 
X-boson momenta into Euclidean space and splitting it into two-parts separated by an 
Euclidean cut-off fi. 

The large momentum part of this integral [between \x and oo] can be simply calculated. 
The large X-boson momentum must flow back via quark- or gluon-lines and can be cal- 
culated perturbatively. The part where the momentum flow-back is through the hadronic 
wave-functions is suppressed by 1/fJ 2 and can be neglected. Precisely as described in 



T5| the result is proportional to cts(fi) and is thus suppressed by 1/N C compared to the 



tree-level matrix elements or it is down by an extra factor of a and reads 



(out | Xj-exchange | in) Short - Distancc = ^ 
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x log— - (o^|X r exchange|m) Leadingl/JVc . (17) 



This part of the integration can thus be included by replacing in Eq. ( |IBJ ) everywhere 
log Mx by log \i leading to 




E 



x , a> s {nc) ( , // . _ . He . 

Oij H Jij log h 7ij log h r {j - r i:j 



+- % log A + % log ^ + r& - rtj 

IT \ J /J, R J (1 J J ) 

+0(a s (ti R )-a s (nc))}C j (ti R ). (18) 

It is obvious from this that the final answer for the decays will not depend on the artificial 
choice of the mass Mx- 

The low energy part of the integral [between and fi] over the X-boson momenta 
must now be performed in a particular low-energy approximation. We must therefore 
identify correctly the vector, axial-vector currents and scalar and pseudoscalar densities 
appearing in Eq. (|13| ) in the low-energy model. This is no problem for the vector and 
axial-vector currents since they are currents from a conserved symmetry. For the scalar 
and pseudoscalar densities the situation is somewhat more uncertain. There is remaining 
dependence on log/ic from the terms with 7^ and 7?-. The scale \ic is the one at which 
currents and densities in QCD and in the low-energy model (CHPT and/or ENJL) are 
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Figure 1: The lowest order, e 2 p°, diagrams for the contributions from Q s . is an 
insertion of the pseudo-scalar currents. The dashed line is an X-boson and the dots are 
the X-boson-pseudoscalar couplings. Full lines are pseudo-scalar propagators. 

matched. We stress again that this problem of matching of the (pseudo) scalar densities 
to the low-energy model, though uncertain, is a much more tractable problem than the 
original matching of the four-quark operators of Eq. (|TTD directly. Practically, we do this 
matching when substituting the CHPT or ENJL quark-condensate by the QCD value at 
He — 1 GeV. For other scales we run the quark density with its QCD anomalous dimension. 



4 Long-Distance 



As described before in [12, [L4J] we can calculate the very low-energy part using CHPT to 



lowest order for the X-boson couplings.^ For the intermediate energy region one can use 
hadronic models for the X-boson couplings. They must have the correct chiral symmetry 
and reproduce at very low energy the model-independent CHPT results. Here we will use 
the ENJL model. As we noted earlier in |TJ|, [L5[ this provides a substantial improvement 
over using lowest-order CHPT X-boson couplings. It brings in no new undetermined pa- 



rameters, describes a fairly wide range of low-energy hadronic phenomena [22] and satisfies 



some short- distance constraints from QCD as well p3 

The contributions to order e°p 2 were already worked out in WA 



We now add the results 
relevant to order e 2 p°. They can be calculated using our earlier methods, a well-defined 
off-shell K-7T Green function or with a reduction method as was used in [ I7| . The 
comparison between both and a more extensive discussion of the results in this section will 
be presented elsewhere [|4| . 

The coupling GeFq has only contributions from Q 7 and Q§. The Wilson coefficients 
C7 is order 1 in 1/N C and Cg is order 1/N C . 

The leading in 1/N C contribution to the matrix element of comes from the two 
diagrams in Figure |l]and is of order N 2 . Using quark- loop diagrams one can prove model- 
independently [p3J] that in the chiral limit the next-to- leading order in 1/N C vanishes to all 
orders in X-boson couplings. The non-factorizable contributions from Qg thus only start 
at order e 2 p 2 and do not contribute to Ge- This contribution leads to GeF$ of order N c . 

4 Notc that we use lowest order CHPT in two different meanings. Lowest order in CHPT in the 
couplings in the Lagrangian of Eq. (ph and lowest order in CHPT in the X-boson couplings which are 



used to calculate the coefficients in Eq. (||). 
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Figure 2: The lowest order, e 2 p°, diagrams for the contributions from Q 7 but next-to- 
leading in 1/N C . Notation is as in Fig. |l|. 

The leading contribution in 1/N C from Q 7 is factorizable but is only of order e 2 p 2 . 
The next-to-leading contribution in 1/N C contributes to order e 2 p° at order N c to the 
matrix-element and leads again to GeFq of order N c . The N c counting presented here is 
scheme-independent. 

We can now calculate the leading contributions using lowest-order CHPT for the X- 
boson couplings and obtain 

SGM = -A lW = CM (19, 

where the Wilson coefficient is the one appropriate for our X-boson scheme as discussed 
in Section ^| and F 2 B (fi) = —(0\qq\0)(fi) the quark condensate in the chiral limit. There 
are also 0(e 2 ) contributions from Q$ to G$ and G27- 
The calculations from the diagrams in Fig. |2| lead to 

e 2 G E [Q 7 ] = -y J^pI C M ■ (20) 

Both results in Eq. ( p!9[ ) and (|20| ) agree with those of ]T^| when the results for 6 2 Ge[Q7] 
listed there are restricted to lowest order CHPT only. Notice that the formula (20) in that 
reference is only valid for A 2 >> M|, My and can only be applied at much higher scales. 
That approximation partially explains the large scale dependence found below 1.3 GeV in 
fT7|j . A there corresponds to /1 here. 

We stress again that the Wilson coefficients in (|T9"|) and (|20|) are in the X-boson scheme 
where the matrix-elements can be identified unambiguously. 

We can now turn to the ENJL model for the X-boson couplings. These results have 
a better high energy behaviour. As mentioned already above, for Q 8 the next-to-leading 
order in 1/N C vanishes and the ENJL result thus coincides with Eq. ([19]). The result for 
e 2 G e[Qt\(ij) I C 7 (ii) is given in Table |I[ The numbers should be compared with the lowest 
order result for e 2 G E [Q$] ENJL = -5040C 8 (/i) using B™ 3L = 2.8 GeV. The latter needs 
to be corrected for the QCD value of Bq(1 GeV) = 1.75 GeV [25 1 and QCD running. This 
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Table 1: The ENJL model results for the contribution from Q 7 to c 2 Ge{^) in the X-boson 
scheme. 

correction is the one discussed in Section ||] in order to identify currents and densities in 
QCD and in the low-energy model and is unambiguous. The usual assumption on the 
value of m s needed shows up in our way in the value of Bq, i.e. the quark condensate in 
the chiral limit. The relevant quantity is really the value of the quark condensate^ rather 
than the value of m H . 



5 Numerical Input 

The input values we use are 

m t {m t ) = (165 ± 5) GeV M w = 80.3945 GeV M z = 91.1872 GeV 
sin 2 (e w ) = 0.2315 m b (m b ) = 4.4 GeV m c (m c ) = 1.23 GeV 

\V US \ = 0.2196 \V cb \ = 0.040 ± 0.002 \V ub /V cb \ = 0.090 ± 0.025 (21) 

Quark masses are in the MS scheme, Mz and Mw are pole masses. The electroweak input 
values have been taken from the review by Sirlin |2"7|] , the CKM matrix values from the 
review by Falk p8| . 

We will simply take sin(5) = 1 so that 

Imr = -6.72 ■ 10" 4 . (22) 

In addition we use two sets of values for the strong coupling constant. Set I we choose 
with the central value of the measurement at the r-mass, as{m T ) = 0.345 [2J| and set II 



with the world average value at Mz, as(Mz) = 0.1186 [27j. We use exact running with the 
two-loop beta function, this is what is most closely needed for the running of the Wilson 
coefficients. 

The scale-independent coefficients 7)1=1,3 needed for the calculation of Ek we take from 
||3(i|l , for fji we interpolate the values given in Table 5 of that reference and obtain 171 = 1.93 
(Set I) and 171 = 1.53 (set II). The others are 772 = 0.57 and 7)3 = 0.47 

In addition we use F K = 112.7 MeV, B K = 0.77±0.10 WB, F = 87 MeV, 5^ CD (1 GeV) 



1.75 GeV and the parameters of the ENJL model as determined by the fit in |E2 



In fact, the next-to-leading order CHPT corrections proportional to quark masses are not the ones 



necessary to change (0|gq|0) into (0|ss|0) 14, g6|. So the matrix element of Qe and Qs at higher orders 
are not proportional to the strange quark condensate or the strange quark mass via PCAC. Interpreting 
the value of Bq in terms of m s via lowest order PCAC as done in the usual £' k /ek analyses it corresponds 
to m s (l GeV) » 140 MeV. 
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6 Results for G$, G27 and Ge 



We now combine the matrix-elements and the Wilson coefficients in the X-boson scheme 
to obtain the quantities needed for the AI = 1/2 rule and e' k /sk- 

We use here the input values given above and the two sets of values for as, labelled I 
and II. In order to have some estimate of higher order corrections in aj we have treated 
the NLO running of the Wilson coefficients and the transition to the X-boson scheme in 
two ways which only differ in higher orders in as- One where we use the exact solution of 
the two-loop evolution equations and multiply with the corrections to go from the NDR to 
the X-boson-scheme. This is labelled with mul for multiplicative in the figures. The other 
option is to only approximately solve the NLO running as a series in as- We also include 
the scheme-correction now additively with the running. This is labelled add in the figures. 
Other estimates give a similar variation. 

For definiteness we also set all the three scales //r, fic and [i equal to each other. The 



variation with this assumption for Bk was shown in |lg] and is smaller than the difference 
between the two values of as used. 

The results relevant for the AI = 1/2 rule are shown in Fig. [| For ReGg we obtain a 
reasonable matching and a value which is in reasonable agreement with the experimental 
value of 6.2. It is dominated by the contributions from Qi and Qi- A large part can be 
considered as the low-energy Penguins as we have discussed in jTJJ. For G 2 7 we do not 
have good stability. The ENJL model improves very strongly on just using CHPT for the 
X-boson couplings as used in the work of the Dortmund group, but it still provides a too 
strong suppression. The curve labelled "I quad" in Fig. ^ shows the CHPT approximation 
for the X-boson couplings. The ENJL improvement is obvious. 

The parts multiplying r of Gg and c 2 Ge are shown in Fig. |j. Notice that the strong 
contribution proceeding mainly through Q e and dominating G 8 is from the gluonic Penguin. 
The electroweak Penguins contribute via c 2 Ge and the contribution is dominated by Qg. 
The stability of both results is quite acceptable. These numbers are the input used for e' K 
below. 



7 B- Parameters in the NDR Scheme 

The matrix-elements are usually quoted in terms of the bag-parameters. These are defined 
as the ratio of the matrix-elements divided by the vacuum- insert ion- approximation (VIA) 
matrix-elements. The latter have no scheme- dependence and so the bag-parameters are 
also scheme-dependent. From our earlier results the bag-parameters can be most easily 
determined in the X-boson scheme. Afterwards they need to be converted to the more 
usual NDR-scheme. 

The usual definitions as given in e.g. || [T7j are 
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0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 



(x [GeV] n [GeV] 

Figure 3: The real part of G$ and G27 as a function of //. The experimental value of 
ReGg is about 6.2 and for G 27 about 0.48. I and II refer to the two values of as used and 
the two curves to two ways of solving the evolution equations as described in the text. In 
the large N c limit both G 2 7 and ReG 8 are equal to 1. 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 

li [GeV] 11 [GeV] 

Figure 4: The part of Gg and c 2 Ge multiplying r as a function of /i. I and II refer to 
the two values of as used and the two curves for each value of as are for the two ways of 
solving the evolution equations as described in the text. 
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(1/2), , 1M [Q 2 



5 A 



(3/2)/ \ n(3/2)/. A M 2 [Ql] 



(J 2)/ , M [Q 6 ] 

c ' 

.12), , Q M [Q 7 ] . 



5f /2 V) = 3v/2- 



I? 

(1/2)/ s M [Qs] 



V 



with 



and 



4 3/2 V) = V2^fi; (23) 



M 2 [Q 3 ] = ((^) I=2 \Q 3 \K°) . (24) 



A = -^(ml-ml: 



C . -16V3L 5 (M P )«M(<-^) 



2? = -2v ^ U|WI 7 w (25) 



(0|gg|0> 2 (//) 

We have restricted the definition of || here to the lowest order in CHPT for consistency. 
Full expressions away from the chiral limit are in the first reference of 0. 
In the large iV c limit one gets the model independent results || 



£(1/2) 


= 3; 




B ( 2 1/2) 


6 

5 ; 




S (3/2) 


= ^ 3/2) 


3 

~ 4 ' 


o(l/2) 
-°6 


_ d(3/2) 
— -°8 


_ R (l/2) 



i ; (26) 

and 

B (3/2) = B (l/2) = Q (2?) 

Again these results are valid in the chiral limit. The operators Q7 and Q$ are discussed 
more below. 
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We normalized the By parameter to the Fierzed part of the vacuum insertion approxi- 
mation. This definition of the B? parameter is the same as used by the lattice community 
3g, and in 



^3|. Notice that in some cases terms beyond lowest order are included in 
their definitions. 

To lowest order in CHPT, i.e. O(e p 2 ) or 0(e 2 p 2 ) for the weak lagrangian one can 

in the notation used here and in 



rewrite the definitions of Eq. 
elements: 



141 for the matrix- 



? (l/2) 



R (3/2) 

-°i* 



(A*) 



5d(/i) 
3 

25C 2 (/x) 
(A*) 



d(3/2) 
-D 2x 



(9G 8 [g 1 ](/i) + G 27 [Q 1 ]( / u)) ; 

(9G 8 [g 2 ](/i) + G 27 [Q 2 ](/i)) ; 
3 



4CM 



3F n 6 



G 8 [Q 6 ](/i) 



"6 X 



80(0|gg|0) 2 ( y u)L 5 (M p ) C 6 (^) 



5 



5 ; 



(3/2) 
7X 

(3/2) 
8x 



(A*) 



5. 



(1/2) 
7X 

(1/2) 
8X 



(A*) 



3F 6 e 2 G E [Q 7 }(,) 



5(0|gg|0) 2 (/i) C 7 (/i) 

e 2 G E [Q 8 ](/i) 



F 6 



5(0|gg|0) 2 (/i) C 8 (/i) 



(2* 



The others can be derived similarly but we have only given the numerically important 



ones and B^ 2 ^ in order to compare with other work 111 



331. We have also restricted 



the VIA matrix-elements to order e°p 2 and 



e 2 p° 



in order to be consistent with our full 
calculation. 

Notice that even in the X-boson scheme the values of B$ , B^\ and B$ depend on 
the value of as because of the running of B (fj,). The bag-parameters B^ and B^ the 
value at leading order in 1/N C and the VIA value, run in the same way with B (fi) [ pij] . 

The large N c result for the Qi operator is 0(e 2 p 2 ) in CHPT and thus vanishes in the 
chiral limit. The lowest e 2 p° order and numerically larger contribution is NLO in 1/N C and 
survives in the chiral limit. This operator together with Q 8 have been studied by several 
groups and techniques and deserves special attention and will studied within the present 



approach in f24 



The numerical results for the bag-parameters in the V-boson scheme can be easily 
derived using the results from Section [| and Ref . |14| , in the case of lowest order CHPT 
couplings of the V-bosons using Eqs. (0) and (p0|) and Eqs. (5.2) and (5.17) of Ref. [|14 |. 
In the case of the ENJL model couplings for the X bosons one needs to use Eq. (|I9"D , Table 



Q and Tables 1, 2 and 3 of Ref. 0. Notice that B. 
at NLO in 1/N C while B { ~ 



(I)X 



1 in the V-boson scheme even 



(i)x 



2.2 ±0.5 26 



In order to get to the NDR-scheme we multiply with the combination — discussed 
in Section H via 



B 



(I)NDR 



1 



\Qi)ln 



E 



3 L 



TT 



ij j 



-(rf. 



TT 



'J 



B? )x (Qj)l 



3 /VIA 



(29) 
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Figure 5: The bag parameters in the NDR scheme for the two values of as we use. 
-^ix nas been divided by 3 to make it fit better in the plot. Notice the large values 

for B^ 2 ^ NDR and B^ 2 ^ NDR responsible for the AI = 1/2 rule and the large value for 



x 

B^ 2 ' NDR which is the reason for our large value of e' k /ek- 



( dx 

The Bj -parameters are in the X-boson scheme. We show the resulting numerical results 
for the inputs given earlier in Fig. [5]. Notice the large values for b[ 1 J 2 ^ ndr and B^ 2 
responsible for the AI =1/2 rule and the large value for B^ 2 ' NDR which is the reason 
for our large value of e' k /ek while B^J 2 ^ NDR shows no such large enhancement. Notice 
though that we get B^J 2 ^ NDR (fi) « 1.3 ± 0.2 for fi between 0.6 GeV and 1 GeV which is 
around 30% to 40% larger than most other analysises. 

In the NDR scheme and in the chiral limit we get B^ 2 [fi) « 0.3 ± 0.2 for scales 
\i between 0.6 GeV and 1. GeV. 



8 Results for ek and e' K to Lowest Order 

The indirect CP- violation as described by ex, defined by Eq. (Q), can be rewritten to first 
order in CP violating parameters, using the AI = 1/2 rule and with Rea >> Ek^^o, 
to a very good approximation as 



e i7r/4 fi m M 12 ^ Ima 



M 



12 



y/2 V Aitlk Reao. 

'(A:) 2 r)!5(x c ) + (X* t )%S(x t ) + 2\* c \* t r) 3 S{x c ,x t ] 



^f 2 K B K m K M 2 w 



(30) 
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Putting^ in the central numerical values of Section |5], sm(8) = 1, ReGj 
0.48 we obtain with c 2 Ge and ImG$ calculated in Section || 



6.2 and G27 



\e K \ = (2.40 - 0.34) ■ 10~ 3 
\e K \ = (2.49 - 0.28) ■ 10" 3 



2.06-10 
2.21 ■ 10" 



-3 



Set I; 
Set II 



(31) 



at a value of /x = 0.75 GeV. The /i-dependence is very mild. The first (second) number 
in brackets is the first (second) term in Eq. (0). The main difference is in the value of 
fji. The effect from the second term which is more dependent on the quantities calculated 
here, is small. 

Given the uncertainty on Bk and the other input parameters this result agrees well 
with the experimental value 



cxp 
£ K 



(2.280 ±0.013) ■ 10" 



(32) 



with $ £K = (43.49 ± 0.08)°. For definiteness we will below use the value of Eq. (|32|) for 
our numerical estimates of e' k /ek- 

Using the same inputs we can now calculate e' K from its definition Eq. (^]). To first 
order in CP violating parameters, using the AI = 1/2 rule and with Reao >> £x ima (b 
we can use to a very good approximation 



1 Rea 2 



Imao Ima2 



With 



K — 1 — 

V2Rea V Rea Rea 2 
tt/2 + 5 2 - 5 = (48 ± 4)° the strong phase 



e k 



(33) 



K 



(9.2 -2.5) • 10" 3 = 6.7- 10" 3 Set I 
(7.4 - 1.9) ■ 10" 3 = 5.5 ■ 10" 3 Set II 



(34) 



at a value of \x = 0.75 GeV. The first (second) number in brackets is the first (second) term 
in Eq. fl33|). The first number is dominated by Q 6 and is usually called the strong Penguin 
contribution. The second one is from the contribution of Qj to Q10 an d is dominated by 
Qs and is usually called the electroweak Penguin contribution. 

We have also shown the /i-dependence of the final result for the two values of cts and 
the two ways of NLO resumming used in Fig. [(]. The predictions are quite stable, partly 
due to the fact that we fully calculate all contributions in the same way. 

Our final prediction to lowest order in CHPT is 



■K 



(6 ±3) • 10" 



(35) 



0(p2) 



6 We use the experimental value of G27 to avoid the large /i-dependence seen in our calculation of that 
quantity. All other quantities were more stable as seen in Section |^. 
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Figure 6: Our result for e' k /ek in the chiral limit using the experimental values of Ek and 
the inputs as described in the text. I and II refer to the two values of as used and the two 
curves to two ways of solving the evolution equations as described in the text. 
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We used the experimental value of ek in the above numbers but using the values we 
calculated would not change any results dramatically. 

The error we quoted there rest on several things, the matrix-elements of Qq and Q& 
depend directly on Bq(1 GeV), this does not affect the relative size of Qq and Q$ but does 
affect the overall value in addition to the typical 30% or so we expect our method to have; 
we did not include any uncertainties in CKM matrix elements and as discussed below, 
there are several more effects which we did not include yet. 



9 Discussion and Corrections 



Our number is significantly higher than any others discussed in the recent literature^]. 
The strong Penguin in our case is enhanced in two ways compared to the usual treatment. 
First the scheme where the transitions to models can be made must be taken into account. 
This provides a first enhancement. Then in the X-boson scheme and in the chiral limit 
the ^-parameter takes a value of B£ x 



2.2 ±0.5 26 



and Bg X 



1 at next-to-leading 

order in l/N c . Putting it together with the first effect leads to _B 6x ndr(/z) « 2.9 ± 0.7 in 
the chiral limit for scales between 0.6 GeV and 1 GeV as discussed above. We do not see 
any similar enhancement for Bg° R (fi) w 1.3 ±0.2, so the usual strong cancellation between 
the strong and electroweak Penguins is much weaker in our calculation. A last reason is 
that we treat both the real and imaginary part at the same order, not as is usually done 
by treating the real part to all orders using experiment and the imaginary part to lowest 
order. 

To the value above one has to add higher order CHPT corrections which are pro- 
portional to pion and Kaon masses. These are of three types, namely: (1) final state 
interactions (FSI) from it — n interactions (2) CHPT corrections which are purely real, 
they are due to AS = 1 higher order terms in the Lagrangian and real loop-diagrams; (3) 
isospin breaking corrections due to quark masses, usually estimated by 7r — r\ mixing, and 
electromagnetic corrections. 

FSI have recently been considered in Ref. |KJ| and later also in [KIT. The main point 



of |34| is that both the real and the imaginary part of the amplitudes should be treated in 
the same way leading to a significant enhancement over the usual estimates. We already 
treat the imaginary and the real parts in Eq. (|33|) in the same way so the main effect of 
FSI is in changing the ratio Rea /Rea 2 from about 16.2 to 22.2. Including only FSI to all 
orders then changes Eq. (|3l|) to 



K 



£k 



(6.8 - 1.8) • 10~ 3 = 5.0 • 10~ 3 Set I 
(5.4 - 1.4) • 10" 3 = 4.0 • 10" 3 Set II 



(36) 



Isospin breaking effects due to the tt° — rj mixing to 0(p 4 ) have been calculated in Ref. 
and references therein. This effect adds a contribution to Im a 2 which is parameterized 
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usually as 



[Im a 2 1 



IB 



n 



IB 



Im ao 



Re ci2 Re ao 

Using the value VLib = 0.16 ± 0.03 |36|] one gets 



(37) 



Ek 



K 



£k 



(6.8 -2.9) • 10~ 3 = 3.9 • 10~ 3 Set I 
(5.4-2.3) • 10~ 3 = 3.1 • 10" 3 Set II 



(38) 



Corrections of type (2) have been fitted to 0{p A 



order in |18| but disentangling this from their FSI effects is not obvious. Notice that in 
and 



in [^] and calculated to the same 

m 



181 the enhancement found for Br is from this source and not from the one discussed 



here. 

Isospin breaking corrections beyond mixing have been discussed recently as well, see 
[38], but their full estimate is not done at the present stage. Our approach allows to include 
them in a more systematic fashion and this will be done in the future. 

Given the uncertainties on the CKM matrix-elements which we did not include and the 
uncertainties on the hadronic matrix elements which we still face, a reasonable guess at 
the final result is 



-K 



(3.4 ± 1.8) • 10" 



This result is somewhat above the present world average |jj 



E K 



cxp 



(2.13 ±0.46) • 10 



-3 



(39) 



(40) 



but quite compatible with it. 

In conclusion we have calculated the chiral limit prediction for e' k /ek in a consistent 
fashion using the 1/N C expansion and the X-boson method. The result is substantially 
larger than most other estimates. After including the two largest expected corrections we 
obtain a result in acceptable agreement with the experiment. 
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A Scheme Dependence 



The matching conditions of Eq. fll5|) between the effective Lagrangian action ( |TTj) and 
the effective X-boson ([T^) in Section [5] include two ten by ten matrices. The explicit 
calculation gives for the gluonic corrections 
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r — r = - x 
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11 
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-11 
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with 



The photonic corrections are 
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(41) 
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(43) 



Ud A - I Ud (AA\ 

n = n u + n d - n = n u -- and n = n u + — (44) 
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and n u the number of up-like light-quarks and rid the number of down-like light-quarks. 
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